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ABSTRACT 

Droplet evaporation is a common phenomenon in 
everyday life. For example, when a droplet of coffee 
or salt solution is dropped onto a surface and the 
droplet dries out, a ring of coffee or salt particles is left 
on the surface(Deegan, et al., 1997). This phenomenon 
exists not only in everyday life, but also in many prac- 
tical industrial processes and scientific research and 
could also be used to assist in DNA sequence analysis, 
if the flow field in the droplet produced by the evapo- 
ration could be understood and predicted in detail. 

In order to measure the fluid flow in a droplet, 
small particles can be suspended into the fluid as trac- 
ers. From the ratio of gravitational force to Brownian 
force a 4 Apg / k B T , we find that particle’s tendency to 
settle is proportional to a 4 (a is particle radius). So, to 
keep the particles from settling, the droplet size should 
be chosen to be in a range 0.1 -1.0pm in experiments. 
For such small particles, the Brownian force will affect 
the motion of the particle preventing accurate meas- 
urement of the flow field. This problem could be 
overcome by using larger particles as tracers to meas- 
ure fluid flow under microgravity since the gravita- 
tional acceleration g is then very small. For larger par- 
ticles, Brownian force would hardly affect the motion 
of the particles. Therefore, an accurate flow field could 
be determined from experiments in microgravity. 

Dcegan et al.(1997) measured the radial, height- 
averaged velocity distribution in an evaporating drop- 
let and observed the formation of a ring. Their experi- 
mental results agreed well with the results of the theo- 
retical analysis which they also carried out. Kan- 
tor(1997) also derived a velocity field using a different 
expression for the evaporation rate than that used by 
Dcegan cl al. Both expressions predict a height- 
averaged outward radial flow in an evaporating drop- 
let. But both theories give only an expression for the 
radial height-averaged velocity. 

In this paper, we will investigate the fluid flow in 
an evaporating droplet under normal gravity, and com- 
pare experiments to theories. Then, we will present our 
ideas about the experimental measurement of fluid 
flow in an evaporating droplet under microgravity. 

THEORY 

Let us consider a droplet on the surface, as shown 
in fig.l. The process of its evaporation is often divided 
into two phases. In the first phase, the contact angle 
decreases while the contact line holds its original posi- 
tion. In an evaporating droplet with a fixed, or pinned, 


contact line, the liquid must flow radially outward for 
the contact line to maintain its position, and the contact 
angle decreases with time as the droplet volume 



decreases. When the contact angle decreases to a criti- 
cal angle, the second phase of evaporation starts. In 
this second phase, the contact line recedes while the 
contact angle remains constant(D. M. Anderson et al., 
1995). The rate of decrease of the contact angle in the 
first phase depends on the evaporation rate of the 
droplet. Therefore, the evaporation of liquid induces 
the fluid flow in a droplet. From the theory of diffu- 
sion and a mass balance in the droplet, the relationship 
between the height-averaged velocity and local evapo 
ration rate was obtained: 


dh _1 d(rvh) 


P H + p r 


dr 


= -J(r,t) 


( 1 ) 


Here, the boundary condition is r=0, rvh=0. v is the 
height averaged velocity, /(r, t) is the local evapora- 
tion rate. Deegan et al. and Kantor respectively devel- 
oped their own expressions for the local evaporation 
rate. From the theory of Deegan et al., the local evapo- 
ration rate is: 

/ = /„(!- rV (2) 


where J 0 is the evaporation rate when X = 0 , and 
X controls the heterogeneity of droplet evaporation. 

From the theory of Kantor, the local evaporation 
rate is: 

(ft • j) = — ^2 arclan^e"^ ^ sinh 2 ( X ) + sin 2 ( jj. ) J (3) 

, x , Z-5283^ 

where A= 17.28(1- " “ 


H)/t exp^ : 


T ) 


H is the 


ambient humidity, S= tan “ ] = tanh(/l ), 0 is the 

\2 / 


contact angle, n is an outward-pointing normal vector 
at some point on the oblate surface, and V and p are 
the parameters of the oblate spheroid, defined such that 
at the surface of the droplet, the r and z coordinates 
satisfy: 
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j + 2 — 1 

(cosh X cos fi) (sinh X sin ju) 

Combining (1),(2) and (3), the radial height- 
averaged velocity v For different theories were de- 
rived. From the theory of Deegan et ah, it is: 


4 t. 


m). 


(4) 


where Af(0), M{t) are the mass of the droplet at 
time 0 and time t, r = is the dimensionless radial 


position of a particle. 

From the theory of Kantor, the radial height- 
averaged velocity is: 



2Ra5 r(1 — r 2 ) 



where a = arctanl 


a 


1 -8 
1 + 5 


and the variable b is 


calculated by integrating the following equation: 


do , A 
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The local radial and vertical velocity V r andv z must 
satisfy the continuity equation: 


dr 


dz 


( 6 ) 


For the solvent with low surface tension, the drop- 
let on the surface is flat, so the lubrication theory can 
be applied. This implies that the radial velocity profile 
is parabolic. We also take the shear stress at the droplet 
ov 

surface to be zero; hence ~~r~ = 0 at z=h(r). Thus, 
az 


v r = A(z-2h)z 


( 7 ) 


The relation between V r and v is therefore: 
_ 3 (lz Z' 

Vr 2 \ h h 'j 

Once the evaporation rate is specified, V r 

are obtained by combining eqns.(4),(5),(6) 
From the theory of Deegan el al., they are: 
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Likewise, from Kantor’ s theory, the local radial 
and vertical velocity V r andV z are: 
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EXPERIMENTS 


( 11 ) 


( 12 ) 


To determine which of the above two theories more 
closely approximates the real behavior, we measure the 
velocity field in an evaporating droplet. In these ex- 
periments, fluorescent particles of diameter around 
0.74p.m are used as tracers to obtain the flow path lines 
in the droplet. 

An inverted fluorescent microscope(Nikon eclipse 
TE 200, 40x Power objective)is used to observe parti- 
cle motion in the droplet. Droplets of radius 800- 
lOOOfim are spotted onto a glass substrate using an 
Eppendorf pipette (adjustable volume:0.5-10pl). 


360 


In order to measure the velocity field in the droplet 
we constructed the simple apparatus shown in fig. 2. It 
is a small cell covered by a millipore membrane to 
slow down the evaporation rate and keep the droplet 
from being disturbed by air currents. The millioprc 
filters used contain 6 x 10 8 holes per cm 2 of diame- 
ter=0.15pm. Using these filters, the droplet evapora- 
tion could be slowed down enough that the droplet 
evaporation time t f could be made as long as 2 to 3 
hours. Such long times are needed to allow time for 
measuring the initial droplet position and size and to 
track carefully a particle position as a function of time. 
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Fig. 3 The surface profiles of the 
droplet at different times 



1. Support 2. Glass coverslip 3. Inverted microscope objective 
4. Droplet 5. Seal Cylinder 6. Porous Membrane 

Fig. 2 Experimental Apparatus 

The suspension of fluorescent particles obtained 
from Polysciences was diluted 2000-fold to obtain a 
solution with a particle concentration of 125ppm. After 
spotting a droplet onto the coverslip and confining it in 
the geometry shown in fig. 2, we first measure the po- 
sition of the contact line along the edge of the droplet. 
We then measure the air-water surface profile of the 
droplet at different times so that the volume vs. time 
can be determined. Finally, we choose a particle in the 
droplet and follow its motion with time, recording the 
particle’s coordinates in the droplet at various times, 
yielding a three dimensional path line. 

RESULTS AND ANALYSIS 

1 . The shape of droplet and droplet evaporation rate 
Since the rate of decrease in droplet volume with 
time is a key parameter in the droplet evaporation 
models, a series of droplet profiles are measured at 
different times. In fig. 3, the symbols are the measured 
droplet surface profiles determined form the maximum 
vertical heights at which particles could be found at 
each radial position r at time t. Data could not be ob- 
tained near the center of the droplet because of the 
limited working distance of the objective, which did 
not permit viewing of the very top of the droplet. Nev- 
ertheless, enough of the profile is obtained to permit a 
fit by the profile expected for the droplet, namely that 
of a spherical cap. 



Fig. 4 The droplet volume as a function of lime 


The droplet volume at each time is calculated using 
the spherical-cap model and plotted in fig. 4. The vol- 
ume of the droplet decreases linearly with lime. This 
implies that the overall evaporation rate of the droplet 
is almost constant, which agrees with the theory of 
Deegan, et al. and the experimental results of Rowan, 
et al.(!995). From theory of Deegan el al., the overall 
evaporation rate rff is: 


_ w)-m 


P 03) 


where V(t) is the volume of the droplet at time t and p 
is the density of the solvent. So, from the experimental 
data of droplet volume vs. time, we can obtain the 
drying time, overall evaporation rate, and initial vol- 
ume of the droplet by fitting experimental data with a 
straight line, as shown in fig.4(solid line). 

The theory of Kanlor predicts an overall evapora- 
tion rate that depends on contact angle as: 


dV 


ol 
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If the initial contact angle is small, eqn. (14) be- 
comes: 
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From eqn. (15), we can see that when the initial con- 
tact angle is small, Kantor’s theory predicts that the 
overall evaporation rate will be almost constant, 

2. Particle’s path in the droplet. 

Two typical experimental results are plotted in figs 
5 and 6. The filled and open symbols are the experi- 
mental radial and vertical positions of a particle as a 
function of time. The radial position(open symbols) 
increases and the vertical position(filled symbols) de- 
creases, showing that the particle moves toward the 
edge of droplet and toward the substrate as drying pro- 
gresses. 
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Fig. 5a The position of particle vs. time for X=0.5 
according to the theory of Deegan, et al. 
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Fig. 5b The position of particle vs. time for H=0.975 
according to the theory of Kan tor 
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Fig. 6a The position of particle vs. time for k=0.5 
according to the theory of Deegan, et al. 
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Fig. 6b The position of particle vs. time for H=0.975 
according to the theory of Kantor 

We also compare the experimental data with the 
theory of Deegan, et al. and of Kantor. Before we cal- 
culate the radial and vertical position of a particle in 
the droplet using the evaporation models, we must to 
find the parameter X in the model of Deegan, et al. 
(equations 9 and 10) and the humidity H in the Kan- 
tor’s model(equations 11 and 12). In order to find X 
and H, we adjust X and H so that the time for a particle 
in the center of the droplet, which is supposed to be 
fixed on the surface of the droplet, to reach the glass 
substrate is equal to the drying time of the droplet. 

Using the fitted X and H, we calculate the radial 
and vertical position of a particle in the droplet pre- 
dicted by the two theories; these predictions are plotted 
in fig. 5 and 6. In each figure the solid line is that cal- 
culated from the theory without Brownian motion, 
while the dashed lines show the range of deviations 
expected to be produced by Brownian motion, using 
the theoretical diffusivity value. In general, the vertical 
position of the particle agrees with the prediction of 
both theories to within the experimental error due to 
Brownian motion. In one case(fig. 5), the predictions 
of the model of Deegan at al. are somewhat better than 
those of Kantor’s model, while the reverse is true in 
fig. 6. The predictions of the radial velocity do not 
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agree as well with the experiments as do the predic- 
tions of the vertical velocity. 

3. Additional predictions 

3.1 The fluid velocity in the droplet. 

From eqns. (9) to (12), we calculate the vertical and 
radial velocities respectively, which are plotted in figs 
7 and 8 for the theory of Deegan et al. and Kan tor at 
times 4400 sec. Here we have taken the droplet radius 
R to be 797pm, its height h(0) to be 365|im and the 
drying time tj to be 4480 sec.. We can see the veloc- 
ity fields for the two theories differ significantly. In the 
theory of Deegan, et al., the radial velocity is higher 
than in the theory of Kan tor, especially near the edge 
of the droplet. Kantor’s theory shows a higher vertical 
velocity than docs the theory of Deegan et al., so that 
when the droplet is nearly dry the velocity vectors 
point almost completely downward; see fig. 8. 



Fig. 7 The velocity field at t=4400s 
according to the theory of Deegan, el al. 



Fig. 8 The velocity field at t=4400s 
according to the theory of Kantor 


3.2 Different evaporation conditions 
By changing X in the theory of Deegan, et al. and 
changing the H humidity in Kantor’s theory, we can 


determine how the evaporation rate affects the calcu- 
lated results. The results are plotted in figs. 9 to 1 2. In 
figs 9 and 10, we can see that the variation in r and z 
becomes larger with increasing X. This implies that the 
more nonuniform the evaporation rate is, the stronger 
the flow is. In figs. 1 1 and 12, a decrease in the hu- 
midity produces an increase in the driving force for 
mass transfer during evaporation, and so the variations 
in z and r also become larger. 



Fig. 9 The effect of X on the time-dependent radial 
position of a particle 



Fig. 10 The effect of X on the time-dependent ver- 
tical(solid lines) position of a particle and the local 
height(dashed lines) of the droplet 



radial position of a particle from Kantor' s theory 
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Fig. 12 The effect of humidity H on the time-dependent 
vertical position of a particle from Kantor’s theory 


CONCLUSION 

1. Experimentally, the overall evaporation rate of an 
evaporating droplet is almost constant. This result 
agrees perfectly with the theory of Deegan, el al., and 
approximately with Kantor’s theory. 

2. The surface shape of droplet is that of a spherical 
cap. 

3. The theory of Deegan et al. and that of Kantor give 
predictions of radial and vertical velocity that arc in 
reasonable agreement with measured velocities, espe- 
cially for the vertical velocity. A precise comparison is 
not yet possible, because of experimental error due to 
Brownian motion. 

4. Near the end of the drying process, the theory of 
Deegan et al. predicts a much higher radial velocity 
than vertical velocity. For Kantor’s theory, the vertical 
velocity is larger than the radial velocity near the end 
of drying. 
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